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ABSTRACT 

We report on correlations between radio luminosity and X-ray timing features in X-ray bi- 
I nary systems containing low magnetic field neutron stars and black holes. The sample of neu- 

tron star systems consists of 4U 1728-34, 4U 1820-34, SerX-1, MXB 1730-335, GX 13H-l,the 
millisecond X-ray pulsars SAX J1808.4-3658 and IGR J00291H-5934, and these are compared 
with the black hole system GX 339-4. The analysis has been done using data from pointed 
f-*) ■ observations of the Rossi X-ray Timing Explorer coordinated with radio observations. In the 

ly^ I neutron star systems the radio luminosity is correlated with the characteristic frequency 

. of the Lfi Lorentzian component detected contemporaneously in the power spectrum, and an- 

ticorrelated with its strength. Similarly, in the black hole system GX 339-4 Lr is correlated 
with the frequency of the Li component in the power spectrum, and anti-correlated with its 
strength. The index of a power-law fit to the correlation is similar in both cases, Lr oc 
O ' and Lr cx {rms)~^^^. At lower timing frequencies, the radio luminosity is further found 

to be correlated with the characteristic (break) frequency of the component of the power 
^2 ' spectra in the neutron stars and, marginally, with the equivalent break frequency in GX 339-4. 
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We briefly discuss the coupling between the innermost regions of the accretion disc and the 
production of the jet and, from the behaviour of the ms accreting X-ray pulsars, the possible 
role of the NS magnetic field. 



binaries: close - stars: binaries : - ISM: jets and outflows radio continuum: stars 



1 INTRODUCTION 

Jet production in X-ray binaries (XRBs) is strongly related to the 
properties of the accretion disc. Recent works have established a 
link between the disc X-ray power and the jet radio power in black 
hole (BH) and neutron star (NS) XRB systems. In BH XRBs a non- 
linear correlation links the radio and the X-ray luminosity when 
in the hard state (e.g. Corbel et al. 2003; Gallo, Fender & Pooley 
2003). Accretion disc and jet theories can translate this coupling in 
to a relation between the power in the jet and mass accretion rate 
(M). Heinz & Sunyaev (2003) have derived a non-linear relation 
between the radio power of the jet observed at a given frequency, 
the mass of the compact object and the mass accretion rate. In the 
case of BHs, it is common use to convert the (bolometric) X-ray lu- 
minosities in Eddington units, and infer the mass accretion rate, that 
is the luminosity is a direct estimator of A/, although, as predicted 
already by e.g. the ADAF model (e.g. Narayan & Yi 1994, 1995; 
Narayan, Garcia & McClintock 1997), observationally the relation 
between luminosity and other AI indicators is not straightforward 
[Homan et al. 2001; see also Homan & Belloni (2005), Remillard 



(2005), McClintock & Remillard (2005)]. In the case of NSs, the 
relation between the X-ray luminosity and other mass accretion rate 
indicators is definitely not one-to-one. Parallel tracks are observed 
in NS XRBs, both for individual sources and across sources, be- 
tween the centroid frequency of the kHz quasi-periodic oscillations 
(QPOs) in the X-ray power spectra and the X-ray luminosity (e.g. 
Mendez et al. 1999; Ford et al. 2000; van der Klis 2001). QPOs are 
thought to be related to accretion disc properties, and in particu- 
lar kHz QPO frequencies are generally interpreted as being related 
to the motion of matter in the accretion disc at a preferential ra- 
dius (see e.g. van der Klis 2004 for a review). This suggests that 
the kHz QPOs may be a more direct indicator, rather than luminos- 
ity, at least of accretion geometry (see van der Klis 2001). Many 
works showed that the characteristic frequency of the kHz QPOs 
are strongly related to the characteristic frequencies of the other 
timing features in the power spectra, so that lower-frequency QPOs 
(much easier to detect than the kHz ones) can be used as a proxy 
(e.g. Psaltis, Belloni & van der Klis 1999; Belloni, Psaltis & van 
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der Klis 2002; van Straaten et al. 2002; van Straaten, van der Klis 
& Mendez 2003; van Straaten, van der Klis & Wijnands 2005). 

A positive correlation between radio and X-ray flux, reminis- 
cent of that found in BH XRBs, has been found in the atoll-type 
NS 4U 1728-34 when in its hard X-ray state [Migliari et al. 2003; 
see Hasinger & van der Klis (1989) for a classification of NSs]. 
Eight atolls and low-luminosity NS XRBs (compared to Z-type 
NSs) have been detected in radio so far, during coordinated radio 
and X-ray observations: 4U 1728-34 when steadily in its hard state 
(Migliari et al. 2003), 4U 1820-34 and Ser X-1 when steadily in 
their soft states (Migliari et al. 2004), Aql X-1 during an X-ray out- 
burst (Rupen et al. 2004), MXB 1730-335 (the Rapid Burster) dur- 
ing X-ray outbursts (Rutledge et al. 1997,1998; Moore et al. 2000), 
the two millisec accreting X-ray pulsars SAX J 1808. 4-365 8 and 
IGR J00291-F5934 during X-ray outbursts (Gaensler, Stappers & 
Getts 1999; Rupen et al. 2002; Pooley 2004). The 'peculiar' atoll 
source GX 13+1, difficult to classify through X-ray analysis, hav- 
ing some X-ray properties of Z-type and some of atoll-type sources 
(but more likely identified as a bright atoll source: Schnerr et al. 
2003), has also been detected at radio wavelengths, and shows a 
radio luminosity similar to (i.e. as strong as) the one of Z sources 
(Roman et al. 2004). 



2 OBSERVATIONS AND DATA ANALYSIS 

We have inspected all the available X-ray observations with the 
Proportional Counter Array (PCA) on-board the Rossi X-ray Tim- 
ing Explorer (RXTE), coordinated with radio observations (see Ta- 
ble of atoll- type and ms accreting X-ray pulsars with a detected 
radio counterpart, i.e. 4U 1728-34, 4U 1820-34, Ser X-1, MXB 
1730-335, GX 13+1, SAX J1808.4-3658 and IGR J00291+5934, 
and of the BH GX 339-4 in the hard state from Corbel et al. (2001). 
Throughout the paper we will use a 'broader' definition of atoll 
sources as (non pulsating) low-magnetic fields NS XRBs accret- 
ing at low mass accretion rates (compared to Z-type NSs), which in 



Table 1. Name of the source, date of the beginning of the observation and 
Obs. ID of the PCA/RXTE observation 



Source 


Date 


Obs ID 


NEUTRON STARS 


4U 1728-34 


2000 May 05 


50023-01-14-00 




2000 May 13 


50023-01-16-00 




2000 May 21 


50023-01-19-00 




2001 May 29 


60029-02-02-00 




2001 lun 01 


60029-02-03-00 




2001 lun 03 


60029-02-04-00 




2001 lun 05 


60029-02-05-00 




2001 lun 07 


60029-02-06-00 




2001 lun 09 


60029-02-07-00 


4U 1820-30 


2002 Jul 25 


70030-03-01-01 


Ser X-1 


2002 May 27 


70027-04-01-00 


MXB 1730-335 


1996 Nov 6 


20093-01-01-00 


SAX Jl 808.4-3658 


1998 Apr 27 


30411-01-10-01 




2002 Oct 16 


70080-01-01-04 




2002 Oct 18 


70080-01-02-01 


IGRJ00291+5934 


2004 Dec 6 


90425-01-01-02 


GX 13+1 


1999 Aug 1 


40022-01-01-01 




1999 Aug 4 


40022-01-02-000 


BLACK HOLE 


GX 339-4 


1997 Feb 4 


20181-01-01-00 




1997 Feb 11 


20181-01-02-00 




1997 Feb 18 


20181-01-03-00 




1999 Mar 3 


40108-01-04-000 




1999 Apr 2 


40108-02-01-00 




1999 Apr 22 


40108-02-02-00 




1999 May 14 


40108-02-03-00 
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this specific case includes also MXB 1730-335 (the Rapid Burster). 
Therefore, hereafter we will devide the sample in atolls and ms ac- 
creting X-ray pulsars. This is a more simplified definition (see van 
der Klis 2005 for a more detailed classification), but it is appro- 
priate for the discussion in this paper. All the observations in our 
sample, with dates, radio flux densities at 8.5 GHz and the best- 
fit values of the characteristic frequencies of X-ray low-frequency 
features are shown in Table|2| 

2.1 X-ray timing analysis 

For NSs, we have used event data with a time resolution of 125 
^s (E_125us_64M_0_ls) for the production of the power spectra of 
4U 1728-34, 4U 1820-30, Ser X-1 and IGR J00291-^5934, and of 
16 ^s (E.16us_64M_0.1s) for MXB 1730-335 (the Rapid Burster) 
and SAX J1808.4-3658. We used time bins such that the Nyquist 
frequency is 4096 Hz. For each observation we created power spec- 
tra from segments of 16s (of 256s for IGR J0029 1-1-5934) length us- 
ing Fast Fourier Transform techniques (van der Klis 1989 and ref- 
erences therein) and we removed data drop-outs and X-ray bursts 
from the data, but no background subtraction was performed. No 
deadtime corrections were done before creating the power spec- 
tra. We averaged the Leahy-normalised power spectra (Leahy et al. 
1983) and subtracted the predicted Poisson noise spectrum apply- 
ing the method of Zhang et al. (1995), shifted in power to match 
the spectrum between 3000 and 4000 Hz. We converted the nor- 
malisation of the power spectra to squared fractional rms (e.g. van 
der Klis 1995). 

For the BH GX 339-4, we have used the same procedures as 
for NS. We used combined binned (B_250us_2A_0_13.Q), single- 
bit (SB.250us_14_17_2s) and event (E^50us_128M_18_8s) data 
with a time resolution of 250 /^s for the production of power spec- 
tra and a Nyquist frequency of 2048 Hz. We created power spectra 
from segments of 128s, averaged them and subtracted the predicted 
Poisson noise spectrum shifted in power to match the spectrum be- 
tween 1000 and 2048 Hz. 

We have fitted all the power spectra with a multi-Lorentzian 
model (e.g. Belloni et al. 2002 and references therein), and plotted 
in the uP^ representation (with Pi, the normalized power and v the 
frequency). In this representation the Lorentzian component attains 
its maximum at the characteristic frequency Vmax- The frequencies 
quoted here are therefore all Vmax values, where Vmax = {vo + 
A^)^''^ {uo is the centroid frequency of the Lorentzian and A the 
half width at half maximum), which are comparable to the central 
frequencies values vq for sharp (quality factor Q=-pp^jfjf > 2, 
with FWHM the Full Width at Half Maximum) features (Belloni et 
al. 2002). 

The features in the power spectra have been identified and 
labeled based on the works of Psaltis et al. (1999), Belloni et al. 
(2002), van Straaten et al. (2002, 2003, 2005), Altamirano et al. 
(2005) submitted, Linares et al. (2005) submitted, Klein- Wolt, van 
Straaten & van der Klis (2005), in prep, (see Klein- Wolt 2004). 



2.2 The sample 

4U 1728-34: Migliari et al. (2003) reported a correlation between 
radio flux densities and characteristic frequency of the X-ray low- 
frequency QPOs. We used the values of the radio flux density at 
8.5 GHz (taken with the Very Large Array: VLA) listed in their 
Table 1. The power spectra were fitted using one broad Lorentzian 
to represent the low frequency noise and the break frequency (Lb or 



Lb2), one or two narrower Lorentzians (Lt and L^) below ^ 50 Hz, 
a broad Lorentzian around 100 Hz (LhHz) and narrow Lorentzians 
to fit the kHz QPOs. A typical power spectrum of 4U 1728-34 is 
shown in Fig. 1. The identification of the features has been made 
using van Straaten et al. (2002). 

4U 1820-30: Migliari et al. (2004) reported a radio detec- 
tion (with the VLA) at 8.5 GHz of the source when steadily in its 
soft (lower banana: LB) X-ray state. Only on 2002 July 25, among 
the seven simultaneous radio/X-ray observations of 4U 1820-30 
taken every 2-3 days, the radio source was significantly detected 
(~ 4(T, although an average over the seven days of radio obser- 
vations gives a 5a level detection). The power spectrum of the 
PCA/RXTE observation on this date is shown in Fig. 1. We needed 
three Lorentzians to fit the power spectrum: a broad Lorentzian to 
fit the very low frequency noise (VLFN), Lb2 at about 9 Hz and 
Lb around 23 Hz. The identification of the features has been made 
using van Straaten et al. (2002, 2003, 2005) and Altamirano et al. 
(2005), submitted. 

Ser X-1: Migliari et al. (2004) reported the discovery (with 
the VLA) of the radio counterpart of Ser X-1 at 8.5 GHz, while the 
source was steadily in its soft (LB) X-ray state. The power spectrum 
of the PCA/RXTE observation simultaneous with radio is shown in 
Fig. 1. We have fitted the power spectrum of Ser X-1 with four 
Lorentzians: a broad Lorentzian for the VLFN, Li,2 at ~ 15 Hz, 
Lb ~ 30 Hz and L^ at 70 Hz. The identification of the features 
has been made using van Straaten et al. (2002, 2003, 2005) and 
Altamirano et al. (2005), submitted. 

MXB 1730-335 (Rapid Burster): Moore et al. (2000) reported 
the discovery of the radio (transient) counterpart of the Rapid 
Burster, detected (with the VLA) during two X-ray outbursts in 
November 1996 and June 1997. Due to the large field of view of the 
PCA, during on-axis pointings of the Rapid Burster the estimated 
flux is contaminated by the presence of the positionally near-by 
bright 4U 1728-34. We have therefore analysed slew data of the ob- 
servation of 1996 November 6, when 4U 1728-34 was not anymore 
in the field of view of the PCA. (The slew observations contempo- 
raneous to the other radio detections have not enough statistics to 
identify the components in the power spectra.) We have fitted the 
power spectrum with four Lorentzians: a broad Lf,2 at ~ 8 Hz, and 
three narrower features at higher frequencies: L(, at ~ 17 Hz, L^ at 
~ 60 Hz and hhHz at ~ 135 Hz (see Fig. 1). The identification of 
the features has been made by us, using van Straaten et al. (2002, 
2003, 2005) and Altamirano et al. (2005), submitted, based on the 
similarities of the power spectra to that of other atoll sources. 

SAX J 1808.4-3658: Gaensler, Stappers & Getts (1999) re- 
ported the discovery of a radio transient emission (with the Aus- 
tralian Telescope Compact Array: ATCA)from the millisec accret- 
ing X-ray pulsar SAX J1808.4-3658 during the decay of an X- 
ray outburst on 1998 April 27. Only an upper limit on the radio 
emission (lower than the detection on April 27) have been found 
at 8.5 GHz on April 26, suggesting a re-flaring. The power spec- 
trum of the quasi-simultaneous PCA/RXTE observation of April 
27 allowed to identify Lb at ~ 0.8 Hz, but due to the bad statis- 
tics above 4 Hz, only a broad not reliably-identifiable feature can 
be fitted around 17 Hz (see also van Straaten, van der Klis & Wi- 
jnands 2005). The PCA/RXTE observation on April 26 allows to 
identify both Lb and L^ in the power spectrum. Rupen et al. (2002) 
reported two radio detections at 8.5 GHz with the VLA, during the 
peak of the X-ray outburst on 2002 October 16 and October 18. In 
the power spectra of the simultaneous PCA/RXTE observations on 
October 16 and 18 we identified five and six Lorentzian compo- 
nents, respectively: L62 (only on October 18), Lb, Llf, Lh, Lhhz 
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Table 2. Name of the source, date of the beginning of the radio observations coordinated with X-rays, flux density at 8.46 GHz, characteristic frequency of 
Lj,2, L;, and (i^break^ a^d for the black hole), the distance to the source and the references. Errors on the radio flux is la, errors on frequencies and 
rms are 90% statistical errors. 



NEUTRON STARS 



Source 


Obs Date 


Fs.s (mJy) 


l'b2 


I'b 




D(kpc) 


Ref. 


4U 1728-34 


2000 May 05 


0.6 ±0.2 


16.14 ±3.82 


22.35 ±0.72 


48.12 ±2.07 


4.6 


M03,GO3 




zuuu jviay vd 


U.oo it U.Uo 


0. 1 U it U.DO 


1 c; Ti _|_ n cT' 


OT Q/1 -1-1 

Z / .04 it 1.41:0 








2000 May 21 


0.62 ±0.10 


7.16 ±0.49 


18.43 ± 0.47 


30.14 ± 1.68 








2001 May 29 


0.11 ±0.02 


— 


2.67 ±0.23 


15.64 ± 0.53 








2001 Jun 01 


0.09 ± 0.02 




2.3d ± 0.58 


13.13 ± 1.45 








2001 Jun 03 


0.11 ±0.02 




1.50 ± 0.23 


13.89 ± 1.19 








2001 Jun 05 


0.15 ±0.02 




1.62 ±0.18 


10.95 ± 0.55 








2001 Jun 07 


0.16 ±0.02 




1.94 ±0.28 


11.78 ± 0.71 








2001 Jun 09 


0.09 ±0.02 




1.37 ±0.11 


8.37 ±0.36 






4U 1820-30 


2002 Jul 25 


0.138 ± 0.035 


9.15 ± 2.38 


23.42 ±0.75 




7.6 


M04,H00 


SerX-1 


2002 May 27 


0.076 ± 0.015 


15.50 ± 1.83 


30.56 ± 0.76 


67.07 ± 2.11 


12.7 


M04,JN04 


MXB 1730-335 


1996 Nov 6 


0.370 ± 0.045 


8.08 ± 4.04 


17.21 ± 1.13 


60.08 ± 2.71 


8.8 


MOO, K03 


SAX Jl 808.4-3658 


1998 Apr 27 


0.8 ±0.18 




0.814 ±0.100 


a 


2.5 


G99,Z01 




2002 Oct. 16 


0.44 ± 0.06 


4.50 ± 1.72 


14.62 ± 0.44 


76.30 ± 2.03 




R02 




2002 Oct. 18 


0.30 ± 0.08 




4.53 ±0.55 


21.43 ± 0.40 






IGRJ00291+5934 


2004 Dec 6*" 


0.180 ± 0.035= 


0.052 ± 0.006 


0.69 ± 0.08 


4.86 ± 0.94 


4 


G05,F04 


GX 13+1 


1999 Aug 1 


0.25 ± 0.05 




4.45 ± 0.66 




7 


GS86,H04 




1999 Aug 4 


0.9-5.2<^ 




11.86 ±1.68 


56.86 ± 2.63 






BLACK HOLE 


Source 


Obs Date 


Fs.s (mJy) 


^break 






D (kpc) 


Ref 


GX 339-4 


1997 Feb 4 


9.10 ±0.10 


0.102 ±0.013 


0.50 ± 0.02 


3.89 ±0.09 


> 7 


C00,H04 




1997 Feb 11 


8.20 ±0.10 


0.107 ±0.014 


0.43 ± 0.02 


3.67 ± 0.07 




Z04 




1997 Feb 18 


8.70 ± 0.20 


0.085 ± 0.010 


0.42 ±0.01 


3.63 ± 0.06 








1999 Mar 3 


5.74 ± 0.06 


0.095 ± 0.009 


0.60 ± 0.01 


4.32 ±0.13 








1999 Apr 2 


5.10 ±0.07 


0.080 ±0.015 


0.35 ±0.02 


3.11 ±0.14 








1999 Apr 22 


3.11 ±0.04 


0.037 ±0.011 


0.21 ±0.06 


2.57 ±0.07 








1999 May 14 


1.44 ± 0.04 




0.05 ±0.01= 


1.41 ±0.12 







a: The power spectrum of this observation allowed to identify Lj, at ^ 0.8 Hz, but due to the bad statistics above 4 Hz only a broad feature can be fitted at 
■~ 17 Hz (see van Straaten, van der Klis & Wijnands 2004). b: The date refers to the start of the radio observation, the public RXTE observation we have 
analysed started on December 7, a few hours after the end of the radio observations, c: The flux density of 0.250 ± 0.035 mJy at 5 GHz (Fender et al. 2004) 
has been converted to a flux density at 8.5 GHz assuming an optically thin spectrum with spectral index a = —0.6 (S^ oc j/", where Sy is the radio flux 
density observed at a frequency v). d: The radio emission during the flare has been converted to 8.5 GHz assuming an optically thin radio spectrum with 
a = —0.6 (e.g. Cooke & Poimian 1991; Penninx et al. 1993; Hjelhning et al. 1990a,b; Fomalont, Geldzahler & Bradshaw 2001). e: the identification is 
ambiguous between L/, (more Ukely) and Inbreak- Ref.: M03=MigUari et al. 2003; G03=GaUoway et al. 2003; HOO=Heasley et al. 2000; M04=MigUari et al. 
2004; JN04=Jonker & Nelemans 2004; MOO=Moore et al. 2000; K03=Kuulkers et al. 2003 and references therem; G99=Gaensler, Stappers & Getts 1999; 
Z01=in 't Zand et al. 2001; R02=Rupen et al. 2002; G05=GaUoway et al. 2005; F04=Fender et al. 2004; C01=Corbel et al. 2000; H04=hynes et al. 2004; 
Z04=Zdziarsky et al. 2004. 



and L„. The identification of tlie features has been made using van 
Straaten et al. (2005) and Linares et al. (2005), submitted. 

ICR J0029 1+5934: this milhsec accreting X-ray pulsar (Eck- 
ert et al. 2004; see Markwardt, Swank & Strohmayer 2004; Gal- 
loway et al. 2005) has been detected in radio, with the Ryle Tele- 
scope, at 15 GHz with a flux density of 1.1 mJy almost at the 
peak of the X-ray outburst (Pooley 2004). The radio emission 
faded rapidly in the following days (Fender et al. 2004), suggest- 
ing an emission from a discrete relativistic outflow, launched likely 
at the peak of the X-ray outburst. We have analysed the power 
spectrum of the PCA/RXTE observation of IGR J0029 1+5934 on 
2004 December 7, the first public observation available with quasi- 
simultaneous radio observations. We could fit the power spectrum 



with three Lorentzians: Lj,2 at ~ 0.05 Hz, at ~ 0.7 and L^, at 
~ 5 Hz (see Fig. 1). The identification of the features has been 
made using van Straaten et al. (2005) and Linares et al. (2005), in 
prep. 

GX 13+1: the 'peculiar' and persistently bright atoll source 
GX 13+1 has been observed simultaneously in X-rays (with the 
PCA/RXTE) and radio at 5 GHz (with the VLA) by Homan et al. 
(2004). They reported a relation between the position on the X-ray 
colour-colour diagram (CD) and the radio emission behaviour (i.e. 
flaring emission in the harder state). The classification of GX 13+1 
as atoll- or Z-type NS is controversial. GX 13+1 shows many X-ray 
properties common to both classes, although more similar to atoll 
sources (Schnerr et al. 2003), but radio properties reminiscent of 
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Figure 2. Top panels: Radio luminosity at 8.5 GHz as a function of tlie cliaracteristic frequency for NSs (a) and of U£ for the BH GX 339-4 (a'); Middle 
panels: Radio luminosity at 8.5 GHz as a function of the rms of L/j for NSs (b) and the rms of for the BH GX 339-4 (b'); Lower panels: Radio luminosity at 
8.5 GHz as a function of ui, for the NSs (c) and lytreak the BH GX 339-4 (c'). NSs are 4U 1728-30 (dots), Ser X-1 (open triangle), the Rapid Burster (open 
diamond), 4U 1820-30 (open circle), and the 'peculiar' atoll GX 13-Hl (asterisks; the dashed vertical line indicates the range in radio luminosity observed 
during the flare, while the timing frequency has been estimated averaging the whole X-ray observation; see Homan et al. 2004). In the BH GX 339-4 we 
marked with dots the observations in hard state before the 1998 X-ray outburst, and with triangles the observations in hard state after the outburst. The upper 
limit on the Vhre-ak frequency in panel c' is the characteristic frequency Vf^ (see Table 1). The slope F of the fitting power-laws (solid lines; GX 13-Hl has been 
excluded from the fit, see §|3j are indicated on the top right comer of each panel (see § 3). 



Z sources (Homan et al. 2004). We have converted the radio flux 
densities from 5 GHz to 8.5 GHz assuming a flat spectrum during 
the steady emission in the soft state, and a = —0.6 (S^ cx v" , 
where is the flux density at a frequency v) during the flare in the 
harder state (consistent with observations of radio flares in NSs: 
e.g. Cooke & Ponman 1991; Penninx et al. 1993; Hjellming et al. 
1990a,b; Fomalont, Geldzahler &. Bradshaw 2001). For the iden- 



tifications of the timing components in the power spectra (Fig. 1), 
we followed the detailed timing analysis in Schnerr et al. (2003). 

GX 339-4: GX 339-4 is one of the two BHs to reveal a correla- 
tion between radio and X-ray flux in hard state over a few orders of 
magnitude (Corbel et al. 2001; the other is V404 Cyg: Gallo et al. 
2003), and the only one to have PCA/RXTE observations, simul- 
taneous with ATCA radio observations, which can trace smoothly 
the X-ray timing behaviour of the source in hard state over an order 
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of magnitude in radio luminosity. In order to compare the relations 
between radio emission and X-ray timing among NSs and BHs, we 
have analysed the simultaneous observations of the BH GX 339- 
4 reported in Corbel et al. (2000; see Table 1), during the hard 
state, and which have an optically thick radio spectrum (i.e. steady 
compact jet). The identification of the features has been made us- 
ing Psaltis et al. (1999), Belloni et al. (2000) and Klein- Wolt et al. 
(2005), in prep. 

Note that previous radio/X-ray timing analysis of other BH 
XRBs shows in the 'peculiar' BH GRS 1925+105 a relation be- 
tween the radio states and the QPO centroid frequency, i.e. the fre- 
quency is lower during the so-called radio plateau state where the 
source has a higher optically thick radio flux (Muno et al. 2001), 
and no evidence for correlations in Cyg X-1, which is almost per- 
sistently in the hard state, but steadily at high X-ray luminosties 
close to the the soft state transition (Pottschmidt et al. 2003). 



3 RESULTS 

We find significant correlations between the radio luminosity 
(at 8.5 GHz and calculated using the distances listed in Table |2} 
and the X-ray timing features and in NS XRBs (Fig. 2a,b,c) 
and Lf in the BH GX 339-4 (Fig. 2a',b',c'). and L;, in the NSs 
and Lf and hbreak in the BH have been chosen because they are 
the components that are almost always present in the power spectra 
thus allowing a comparison among observations and sources. 

Spearman rank tests give positive correlations for Lh — Uh 
(99.2%) and Lij — ut (99.1%) without considering ms accret- 
ing X-ray pulsars (i.e. with 4U 1728-34, 4U 1820-34, Ser X-1, 
MXB 1730-335, GX 13+1). Adding also the ms accreting X-ray 
pulsars (IGR J00291+5934 and SAX J1808.4-3658, not shown in 
Fig. 2 and discussed in more detail in S I4.2> we still found a signifi- 
cant correlation for L_r — Vb (99.3%), but only marginal for L_h — Vh 
(97.4%). 

As expected from the negative correlations found between the 
characteristic frequency and the rms of the Lorentzian components 
in the power spectra of NS XRBs (e.g. van Straaten et al. 2002, 
2003, 2005), we also found a negative correlation between and 
the rms of hh (99. 1 % with atoll sources and 99.7% with also ms X- 
ray pulsars) and the rms of L;, (98.8% with atoll sources and 99.7% 
with also ms X-ray pulsars). 

We fitted all these correlations with power-laws. The fits with 
a power-law model give /d.o.f . not better than ~ 4 (considering 
errors on both variables), however, due to a lack of physical models 
to describe such correlations, the derived fitting parameters, like 
the power-law spectral index, are useful to quantify the relations 
of these observational quantities. To quantify the likelihood of the 
power-law model to describe the correlations, we quote also the 
value of the linear coefficient r calculated for each correlation in 
a logarithmic scale. In order to be more conservative, we chose to 
perform the fit considering the atoll sources without the 'peculiar' 
atoll source GX 13+1. In this way we can trace the behaviour of 
the other sources (GX 13+1, but also of the ms X-ray pulsars) with 
respect to them (Spearman rank tests without GX 13+1 and the 
ms accreting X-ray pulsars gives: for L_b — a significance of 
98.7%, for Ljj — a significance of 98.9% and for L_r— rms of 
Lji a significance of 97.7%). In the following we quote the fitting- 
parameters' values with la errors derived using only errors on L_r 
(similar values are derived using errors on both variables). 

Fits with a power-law L_r = y4 x i^^y^ of the radio luminosity 
as a function of the timing frequencies give, for a slope of F = 



1.30 ± 0.10 and a normalization A — 10.8 ± 3.8, and for i/t a 
r = 0.63 ± 0.04 and A = 197.1 ± 17.0. Note that the Ljj - i^t 
correlation looks more 'bimodal' than L_r — (Fig. 2,c; r = 0.85 
and r — 0.93, respectively, for L_r — vt and L_r — Uh). 

A fit with a power-law of Lh with the rms of the timing com- 
ponents give, for the rms of Lj,: F = —2.32 ± 0.14 and A = 
7.303(±2.255) x 10* with a correlation coefficient r = -0.87; 
and for the rms of U: F = -1.66 ± 0.20 and A = 7.65 ± 3.81 
with a correlation coefficient r = —0.66. 

Correlations similar to those found for i^h and for the break 
frequency in NSs hold also for the BH GX 339-4 in the hard 
state. In Fig. 2a' we show the radio luminosity at 8.5 GHz as a 
function of the characteristic frequency in the BH GX 339- 
4 (see § 14.31 for a discussion). The fit with a power-law gives 
A = (4.88 ± 0.11) X 10^ and a slope of F = 1.37 ± 0.02 (slope 
consistent with being the same as with i^h in atoll NSs). The rms of 
he, as in NSs, decreases with radio luminosity and a power-law fit 
gives A = 3.54(±0.45) x lO'' and a slope of F = -2.38 ± 0.03 
(with a correlation coefficient r = —0.90; Fig. 2b'). This slope is 
also consistent with the F found for the rms of Lh in atoll NSs. A 
marginally significant correlation (95%) has been found between 
the radio luminosity and the break frequency vtreak (Fig. 2,c'). A 
power-law fit (without the upper limit point) gives F = 0.87± 0.02 
and A = 2.59(±0.99) x 10^, with a linear correlation coefficient 
of r- = 0.88. 

As a caveat we note that, in NSs, in the case of a transient 
outburst. X-ray and radio emission may have been decoupled. X- 
ray timing features, being related to the disc properties, should be a 
more reliable tracer of the compact steady continuous replenished 
jet, rather than of the evolution of a transient optically thin jet, at 
least after the jet has been launched and thus decoupled from the 
disc. In this context, this caveat specifically applies to the ms accret- 
ing X-ray pulsar SAX J1808.4-3658 detected in radio only during 
X-ray outbursts that do not, in fact, follow the Ln — i^h correlation 
found for the other sources (Fig.|3l see §|42}- The fact that the data 
of the (likely) transient jet in the Rapid Burster (and actually also 
the ms accreting X-ray pulsar IGR J0029 1+5934 which is in the last 
part of the X-ray outburst decay) instead do lie on the radio-X-ray 
timing correlations may suggest a continuity in power between the 
compact and the transient jet at launch (see discussion in Fender, 
Belloni & Gallo 2004 for BHs). 



4 DISCUSSION 

We have analysed all the available PCA/RXTE observations co- 
ordinated with radio, of atoll NSs and ms accreting X-ray pulsars 
with a detected radio counterpart. 

We find significant correlations between the radio luminosity 
at 8.5 GHz and X-ray timing features. We compared these rela- 
tions in NSs with those of the BH GX 339-4 in the hard state, for 
which we could smoothly follow the radio behaviour and the high- 
resolution X-ray timing variations over an order of magnitude in 
radio luminosity. We found that similar relations hold in the two 
different classes of sources. In particular we observe: in atoll NSs 
correlations between radio luminosity and and between radio 
luminosity and the break frequency component L;,, in BHs corre- 
lations between radio luminosity and L< and (marginally) between 
radio luminosity and the break frequency Ltreafc- In the following 
we will discuss some implications of these results and suggest pos- 
sible scenarios, which require testing using a larger sample. 
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4.1 Atoll-type NSs 

The phenomenon of the parallel tracks in the X-ray luminosity vs. 
kHz QPO frequency plot for NSs (e.g. Mendez et al. 1999; see van 
der Klis 2005 for a review) may be explained using a single time 
dependent physical quantity, usually inferred as the mass accretion 
rate in the disc Md, if we consider the kHz QPO frequency to be 
governed by the balance between this quantity and luminosity Lx , 
while the luminosity responds to both this quantity and its time- 
averaged variations (van der Klis 2001). The QPOs seem to be a 
good tracer of the disc geometry as determined by the balance be- 
tween Md and Lx, in the regions very close to the compact object 
where also jets are thought to originate (see also below), van der 
Klis (2001) also suggested that one possible mechanism for bal- 
ancing the energy budget is by removing material from the binary 
in the form of a mildly-relativistic jet, so that the accretion rate on 
to the compact object is no longer the same as the mass transfer rate 
through the disk. 

All the variability components in the power spectra follow a 
universal scheme, when plotted against the upper-kHz QPO. There- 
fore the disc geometry may be inferred also by low-frequency tim- 
ing features. In particular in atoll sources, a tight correlation exists 
between the centroid frequency of the upper-kHz QPO and that of 
(the radio power tracer) L^: the best-fit power-law is Vh,o oc v^ q 
(van Straaten, van der Klis & Mendez 2003), or i/h oc u'^''^ using 
I'max (van Straaten, van der Klis & Wijnands 2005). If we assume 
that Vu is related to an orbit in the disc at an inner radius Ri„, the 
fact that the radio jet power increases with Vh, may be explained by 
a scenario in which the jet particles and magnetic field lines may 
be 'squeezed' as the disc moves inwards. (The jet power possibly 
increases until the magnetic field lines of the jet interact with the 
magnetic field of the NS surface; see also below in i) l4.2t . 



4.2 Millisec accreting X-ray pulsars 

The observations of the ms accreting X-ray pulsar 
IGR J00291-I-5934, taken in the last part of the decay of the 
X-ray outburst on 2004, also lies on the radio luminosity/X-ray 
timing frequency correlations of atoll sources (adding this source 
the Lfl — i^b and La — rank correlations are significant to the 
99.6% and 99.7% level, respectively). On contrary, the observa- 
tions of SAX J1808.4-3658 seem to diverge from the radio/X-ray 
coupling found in atoll sources and IGR J0029 1-1-5934, especially 
from the Lr — Vh relation (see § |3] note that the ms X-ray 
pulsars are still consistent with the rank correlations between 
radio luminosity and the rms of the X-ray timing features of atolls 
sources; see also Fig. 3, middle panel). The radio upper limit of 
SAX 11808.4-3658 has been derived on 1998 April 26 from an 
observation during the last part of the decay of an X-ray outburst, 
as for IGR J00291-I-5934, and this value is also consistent with the 
correlations of the atolls. The observations that do not lie on the 
radio luminosity/X-ray timing frequency correlation are those of 
SAX J1808.4-3658 taken during the peak of the outburst in 2002 
and the one that apparently re-flared on 1998 April 27. Although 
the lack of observations does not allow a complete understanding 
yet of the behaviour of the simultaneous radio/X-ray emission of 
ms X-ray pulsars throughout the evolution of an X-ray outburst, 
based on the present data we can start to make some remarks and 
speculations. 

Focusing on the radio luminosity/X-ray timing frequency re- 
lations (Fig.|3] upper and lower panels) we see that when the fre- 
quency is high {uh > 20; > 4) the radio luminosity is below 
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Figure 3. Radio luminosity at 8.5 GHz as a function of (upper panel), 
the rms of (middle panel) and (lower panel) for the ms accreting X- 
ray pulsars IGR J00291+5934 (open cross) and SAX J1808.4-3658 (stai's). 
The gray region represent the power-law fit to the atoll sources considering 
1(T eiTors on the normalization A only (e.g. 83.5=^ x (v)^ \ see §|3j. The 
question mark indicates an uncertainty on Vy^, because since no Ljj could 
be identified on 1998 April 27, we used the value of the X-ray observation 
closest to that day and for which v\y was consistent with the one on 1998 
April 27 (i.e. group number 21 in van Straaten et al. 2005; see also !; 12.21 . 
The arrows indicate the chronological sequence of the observations during 
each X-ray outburst. 
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the one expected from atoll sources, while it is above during the 
re-flaring when the source show lower timing frequencies. Further- 
more, we see that in the 2002 outburst, in two days the timing fre- 
quencies change of a factor of four while only a little decrease is 
observed in radio luminosity, i.e. an increase in radio loudness. This 
may be an indication of an actual decoupling of disc and jet during 
the outburst, thus signature that a discrete jet has been launched. 
The fact that when the timing frequencies are higher the source is 
less radio loud than at lower timing frequencies, may be due to 
the relative importance of the magnetic field of the NS when the 
disc is getting closer to the NS during the outburst. For the two 
observations in 2002 we have a direct measurement of the upper 
kHz QPOs: L„ ~ 700 Hz on October 16 and L„ ~ 400 Hz on 
October 18 (see also van Straaten et al. 2005). If we relate this 
frequency to the Keplerian orbit in the disc at a radius Ri„, for a 
1.4 Mq NS we derive Kin ~ 70 km in the October 16 observa- 
tion and Kin ~ 100 km in the October 18 observation. For the 
observation of 1998 April 27, using the relation Vh oc v'^^ we de- 
rive Kin ~ 140 km. In ms accreting X-ray pulsars the magnetic 
field of the NS might be higher than in atoll NSs (e.g. Chakrabarty 
2005), and for SAX J1808.4--3658 it has been estimated to be 
~ 10* - lO'' G (Psaltis & Chakrabarty 1999). The (yet poor) ob- 
servations may suggest that below, say, ~ 10*^ km from the NS, 
the magnetic field, interacting more strongly with the regions of 
the disc where the jet is formed - closer to the compact object as 
the jet gets more powerful as suggested by the correlation with the 
X-ray timing frequencies - suppresses in some way the power of 
the jet (possibly by interfering with a magnetic expulsion mecha- 
nism; see also i; l4.U . We can have a rough estimate of the radius 
at which this happens if we assume that this interaction is strong 
when this region is close to the Alfven radius Ka- Assuming for 
SAX J1808.4-3658 a magnetic field at the surface of ~ 5 x 10* G 
(Psaltis & Chakrabarty 1999), a luminosity during the decay of the 
1998 outburst on April 27 of ~ 3 x 10^^ erg/s [from Gilvanov et 
al. (1998), using a distance of 2.5 kpc (in 't Zand et al. 2001)], and 
a radius of the NS of ~ 10 km (Burderi & King 1998), we obtain 
Ra ~ 110 km, and slightly smaller radii for the other observations 
(see e.g. Frank, King & Raine 2002 and references therein). This ra- 
dius is consistent with the inner radius we have estimated with kHz 
QPOs' frequencies, therefore suggesting that the magnetic field of 
the NS, through its interaction with the innermost regions of the 
accretion disc, might play a role also in (e.g. inhibiting) the pro- 
duction of the jet. 

4.3 NSs and BHs 

Casella, Belloni & Stella (2005) found a clear association between 
a low-frequency QPO, the so-called type-C QPO in BHs (e.g. Wi- 
jnands et al. 1999; Remillard et al. 2002; Casella et al. 2004) and 
the horizontal-branch oscillation (HBO) in Z-type NSs, suggesting 
a similar physical origin for these two timing features, van Straaten 
et al. (2003) already compared the low-frequency features of atolls 
with those of Z sources, and, based upon the correlation with the 
frequency of the kHz QPOs, identified the HBO in Z sources with 
Lft in atoll sources. [Note that, a correlation between radio power 
and the position in the CD has been found in Z sources, and it is 
precisely in the HB that Z sources are more radio loud (Penninx 
et al. 1988; Hjellming et al 1990a,b).] In BHs, Psaltis et al. (1999; 
see also Fig. 12 in Belloni et al. 2002) found a tight almost linear 
correlation between the frequencies of the narrow low-frequency 
QPO Llf (i.e. also the type-C QPO) and a broader component that 
they called (reminiscent of the lower kHz component in NSs). 



Using these known relations we can (linearly) associate Lh of NSs 
with Lf of BHs (where the choice for these two features have been 
dictated by availability in the power spectra of our sample; see §|3}- 
The fact that we find for these two features the same scaling with 
radio luminosity supports the idea that the same physics lie behind 
their origin. 

Since all the compact jet models predict that the total jet power 
Lj and the power radiated in radio L_r are related as L_r ocLj* 
(e.g. Blandford & Konigl 1979; Falcke & Biermann 1996; Fal- 
cke, Markoff & Fender 2001), the finding that the radio luminosity 
scales with Uh and ui, in NSs and BHs respectively, with a slope 
of F ~ 1.3 — 1.4, would imply an almost linear correlation of Lj 
with these timing features. 

It is interesting to note that if we consider the frequencies of 
the timing components correlated with the radio luminosity (i.e. Uh 
for NSs and vi for the BH GX 339-4) as originating from Keple- 
rian motion of matter in the disc, from the dynamical timescales 
observed we derive a distance to the compact object of ~ 100 — 
300 Kschvj and 100-500 Kschw (where Kschvj = 2GM/c^ is the 
Schwarzschild radius), respectively in NSs and in GX 339-4. These 
radii are consistent with those derived from high-resolution radio 
imaging of the active galaxy M87 (Junor, Biretta & Livio 1999) 
and consistent with the idea that jets are formed very close to the 
compact object. 

The relations between the radio luminosity and the break fre- 
quency are even more intriguing, because such features have been 
observed also in AGN (e.g. McHardy et al. 2004). Therefore, these 
correlations, if further confirmed especially with a larger sample of 
BHs, open the possibility of the existence of a new 'fundamental 
plane' for XRBs and AGN (Merloni, Heinz & Di Matteo 2003), 
where the dimensions are the mass of the compact object, the ra- 
dio luminosity and the characteristic frequency of a timing feature 
(e.g. break), the last being, contrary to the X-ray luminosity (see 
Merloni et al. 2003 and Falcke, Kording, Markoff 2004), indepen- 
dent from the distance to the source. X-ray timing features can be 
the key features to finally find the common physical link between 
the accretion disc and the jet radio emission in NSs and BHs of all 
masses. 
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